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Abstract. We study the relation between the mean elfective surface temperature Ts and the internal temperature T\, for magnetic 
neutron stars, assuming that the magnetic field near the surface has a presumably small-scale structure. The heavy-element (iron) 
and light-element (accreted) heat-blanketing envelopes are considered, and the results are compared with the case of a dipole 
magnetic field. We argue that the difference in the rb(rs)-relation for different magnetic configurations is always much smaller 
than a possible difference caused by variations of the chemical composition in the envelope. 
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1. Introduction 

The most direct evidence for the nature of a poorly known 
superdense matter in neutron stars is likely to emerge from 
comparisons of surface thermal radiation meas urements with 
I predictions of neutron star cooling models ( e.g.. |Pethickll992t 
" |Page'l997'.'l998':' Yakovlev & Pethick'20o7). ModeUng of neu- 
tron star cooling is a complex problem that, generally, re- 
quires calculations of the temperature profile from the surface 
to the core at various stages of the evolution. The tempera- 
. ture profile for nonmagnetic neutr on stars has been the sub- 
' iect of study for many au t hors ("e.g-loudmundsson et alJl983t 
iNomoto & Ts uruta '1987'; 'Potekhin et alJI 19971 and references 
therein). The key issue of these studies is the so-called 7^(7^)- 
relation used in cooling models (T, and are the surface and 
internal temperature, respectively). 

Most of the neutron stars, however, possess surface mag- 
netic fields B ~ lO'^-lO'^G, and some neutron stars are possi- 
bly magnetars with B > 10^* G. The internal magnetic field can 
be even higher. Such strong fields can affect the properties of 
plasma in neutron stars and alter the rbCr,) -relation. In general, 
the magnetic field strength and direction varies over the stel- 
lar surface, and hence may be different for different surface 
points. In many applications (e.g., in the neutron-star cooling 
theory) it is sufficient to know the mean effective surface tem- 
perature Ts instead of the position-dependent T^. The effective 
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temperature is defined by the Stefan law, L - A-nR crT^ , where 
cr is the Stefan-Boltzmann constant, and L is the thermal lumi- 
nosity in a local neutron-star reference frame, integrated over 
the surface. The apparent luminosity measured by a distant ob- 
server IS, L°° - {I - rg/R)L, and the apparent surface temper- 
ature inferred by the obs erver from the radiation spectrum is 

= r, VI - ^s/R re.iT..lThornell977h. 

The effects of a strong magnetic field on thermodynamic 
and kinetic properties of the outer neutron star layers have been 
reviewed, for instarice, by Yakovlev & Kaminket ( 1994) and 
IVentura & PotekhinI ( 120011) . The thermal structure of magne- 
tized neutron stars has been analyzed by a number of authors, 
often adopting a simplified magnetic geometry. For instance, 
much attention has been paid to the cases of the radial mag- 
netic field (e.g..'Hernaui stll985H Van Ripei*1988^.'Schaaf'l99d 
Hevl & Hernauist 1998) or the tang ential field (with the field 
fines parallel to the surface; e. g.. lHernquistl 1 98.5HSchaai 1 990l: 
iHevl & Hernquistii.l998>) . The quantizing radial field field de- 
creases the difference between and T^, whereas the tangen- 
tial field increases this difference. 

The case of an arbitrar y inclinat ion of the field lines 
has be en considered by Greenstein & Hartke ( 19831) and lPagel 
ill 99.51) who argued that the rb(rs)-relation for such a field can 
be constructed in a simple way from the relations for the radial 
and tangential m agnetic fields. Numerical calculations of the 
therrn al structure jPotekhin & Yakovievll200lUPotekhin et all 
12003') confirm an accuracy of this approximation. These cal- 
culations show that the rb(7^s)-relation for the dipole mag- 
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netic configuration is almost independent of the field for B < 
a few X lO''' G. The reason of such behaviour is the compensa- 
tion of a decrease in the thermal flux near the magnetic equator 
by an increase in the flux near the pole. 

This result can essentially simplify the cooling calculations 
for magnetized neutron stars with dipole magnetic fields. It is 
likely, however, that the surface magnet ic field o f neutron stars 
depart s fr om a dipole configuration. lArons & Scharlemanrj 
ill 9791) and lAronj 1 (199? .' 2000) noted that pulsars with long pe- 
riods require a more complex field configuration than a dipole 
if pair creation is es sent ial for the mechanism of radio emission. 
iGil & Mitra'(2001) and lOil & Melikid ze (2002) also noted that 
radio emission of many pulsars can be explained if one adopts 
the model with a strong (> 10^^ G) and complex surface 
field with a small c urvature of the field lines (~ 10^ cm). 
lGil&Sendv]3 l l2OO0l) found that the behaviour of drifting sub- 
pulses observed in many pulsars is consistent with the vacuum 
gap maintained by a strong sunspot-like magnetic field. Recent 
observations of the X-ray spectra of some pulsars provide an 
opportunity to estimate the magnetic field near the neutron star 
surface. For example, a possible interpretation of the feature 
observed in the spectrum o f PSR B 182 1 -24 as a cyclotron 
emission jBecker et alJ2003[ see, however. iMineo et alJl2004() 
indicates that the local magnetic fields on the neutron star sur- 
face can exceed the conventional dipole field inferred from the 
spin-down data. 

The growing evidence of the distinction between the lo- 
cal field strength at the stellar surface and the global dipole 
field suggests that this can be a general phenomenon in neu- 
tron stars. Several theoretical explanations of this phenomenon 
were suggested. Rudermani til99L) considered a "plate tec- 
tonics" model, in which a complex configuration of the sur- 
face magnetic field of a neutron star is a re sult of neutron- 
star cr ust cracking and platelet movements. iRuderman et alJ 
( Il998l) f ound that this model a grees with observations of pulsar 
glitches. Geppert et al. ( 2003) suggested a Hall-driven mecha- 
nism of the formation of the sunspot-like magnetic-field struc- 
ture at the neutron-star surface. The presence of small scale 
field components can also be generated during the initial con- 
vective s tage of evolution, e.g. , by the turbulent mean-field 
dynam o jThompson & Duncanl Il993t iBonanno et all l2003l 
boos *). The unstable stage in proto-neutron stars lasts ~ 30^0 
s, which is sufficient for dynamo to reach a saturation level. The 
magnetic field generated in proto-neutron stars will be frozen 
into the crust that starts to form almost immediately after con- 
vection stops. Since the crustal conductivity is high, both the 
large and relatively small scale (~ 10^ cm) fields can survive 
during a long time comparable to the lifetime of pulsars, ~ 10^ 
yr cUrpin & Gil 2004). 

If the dipole field in neutron stars is accompanied by a 
stronger small-scale field, then the thermal structure of the sur- 
face layers can be changed qualitatively. For instance, the dif- 
ference between the equator and polar temperature should be 
much reduced. The small-scale field can also affect the total 
photon flux at given T\, - that is the T\,(Ts) relation. In this pa- 
per, we consider the thermal structure and calculate the T\,(Ts) 
relation in the case of a neutron star with a strong small-scale 
magnetic field in the envelope. 



2. Statement of the problem 

The thermal structure of the neutron star envelope is consid- 
ered in steady state plane-parallel approximation. We assume 
that Tb is equal to the temperature at the inner boundary of the 
envelope (taken at the neutron drip density 4x10" gcm"^) and 
does not vary over this boundary. 

The statement of the problem for calculations of the tem- 
perature profile in neutron star envelopes has been described in 
detail by Potekhin et al. (2003). In the present study, the only 
difference is the expression for the heat flux, that should incor- 
porate the effect of small-scale magnetic fields. In the magnetic 
field B, the heat flux is related to the temperature gradient by 

F = -^||(5)V||r - kAbw^t - KAB)b X vr, (i) 

where the tensor and vector components along and across the 
field are marked by || and ±, respectively, and /Ca is the so- 
called Hall component; b - B/B. Generally, all components of 
the thermal conductivity tensor depend on the field strength, B. 
The vector components in Eq. ([0 are 

ViiT ^b-(b-VT), V_lT = VT-b (b VT). (2) 

We assume that small-scale fields are stronger than the 
large-sc ale field (as it follows from the t urbulent dynamo mod- 
els, e.g.. lThompson & Duncantl993;.Bonanno et alJ2005l) . and 
that the length-scale of a temperature variation over the sur- 
face is much larger than the vertical length-scale. This allows 
us to apply the local plane-parallel approximation for studying 
the m ean thermal sti'ucture of the heat-blanketing en velope (see 
iGudmundsson et al]ll983 ' and Pot ekhin et al.l2003[ for discus- 
sion). Then, averaging Eq. ([0 over the directions of the mag- 
netic field (that is approximately equivalent to averaging over 
the stellar surface), we obtain for the mean heat flux 

{F}=1[^(B) + 2kAB)]^, (3) 

where z is the local proper depth iri the env elope. This is equiv- 
alent to Eq. (7) of IPotekhin et all (l2003l) . but wifli flie effec- 
tive thermal conductivity k replaced by its average over all an- 
gles. Thus, unlike Potekhin et al. (.2003) . where the variations 
in strength and direction of the magnetic field were assumed 
smooth over the stellar surface and small compared to the tem- 
perature gradient in the heat-blanketing envelope, which al- 
lowed us to treat the field as locally constant, here we consider 
the opposite limit of a rapidly varying field direction (i.e., "en- 
tangled" field lines). For simplicity, we assume in what follows 
that the average strength of small scale magnetic fields does not 
depend on the depth. 

We calculate the thermal structure of a neutron star en- 
velope by integrating Eq. dSl l with the use of the nume rical 
scheme developed previously ('Potekhin et al.lll997ll2003l) . We 
examine models of envelopes composed of iron and accreted 
material. Th e accreted envelope is assumed to be the same as 
described bv lPotekhin et al.l d 19971) : the surface H layer is fol- 
lowed by He, C, O, and Fe layers. The physics input (equation 
of state, radiative o pacities, and el e ctron thermal conductivi- 
ties) is described in Potekhin et alJ ( l2003l) . Note some uncer- 
tainty of this input for the outermost layer of the iron envelope 
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Fig. 1. The redshifted effective surface temperature versus the internal temperature for a canonical neutron star (M = 1.4 M©, 
R = 10 km) with the magnetic field B = 10'^, 10'^, and 10'"^ G (left, middle, and right panels) and for the case of dipole (dot-dash 
lines) and small-scale (solid lines) fields. The upper and lower curves represent the accreted and iron heat-blanketing envelope, 
respectively. 



where Fe plasma is only partly ionized. In this case, the equa- 
tion of state and opacities are based on the Thomas-Fermi and 
mean-ion approximations. For the accreted envelope, in con- 
trast, we use the accurate equation of state and opacities for the 
partially ionized hydrogen plasma in a strong magnetic field 
jPotekhin & Chabrieii2004j and references therein). 



3. Results 

The calculated temperature profiles and their dependences on 
B and T\, are quite similar to those presented in Potekhin et al. 

Therefore we will not describe the results in detail but 
focus on the comparison of the Tb{Tf;) relation for the small- 
scale and dipole fields. 

The examples of such a comparison are shown in Figs. \l\ 
and 121 The panels (from left to right) in Fig. ^ correspond to 
B = 10** G, 10'2 G, and 10"* G, typical for millisecond pulsars, 
ordinary radio pulsars, and anomalous X-ray pulsars, respec- 
tively; B is the average field strength in the case of a small-scale 
field and the field at the magnetic pole in the dipole field model. 
For convenience, the non-magnetic rb(rs)-relations are shown 
in the middle and right panels of Fig.[nby the dotted lines. 

Figure |2| shows the dependence of on B for different 
models and for two values. Independent of the field geome- 
try, a weak field (B < 10^ G) does not affect the rb(rs)-relation 
for both the iron and accreted envelopes. With increasing B, the 
effective temperature first decreases and then increases non- 
monotonic ally. This behaviour is explained as follows. The 
value of Ts at a given Th is controlled mainly by the opacity 



values in the "sensitivity strip" jGudmundsson et allll983l) lo- 
cated typically at p ~ 10'*-10^ g cm~^, where heat transport 
is provided by the electron conduction. In a strong magnetic 
field, the electron thermal conductivity as a function of p un- 
dergoes oscillations when the electron Fermi energy crosses 
the magnetic Landau levels (see, e.g., 'Yakovlev & Kaminkeil 
11994: Ventura & Potekhin 2QQ1J- The first oscillation at p = 
pB = 7045 Z?j2^ (A/Z) g cm"^ is the strongest one (here B12 = 
B/IO'^ G and A and Z are the atomic weight and charge num- 
bers); the classical (non-oscillating) magnetic electron conduc- 

— 00 

a decrease of first 
occurs in the classical regime, where the effect of the magnetic 
field amounts to a suppression of /c^ in Eq- Q. With further 
increase of B, the field becomes quantizing in the sensitivity 
strip, leading to an increase of . The increasing pieces of 
curves in Fig. |2] are wavy because of quantum oscillations of 
K\\, which are smoothed (but not entirely) by the integration of 
Eq. (|3} over the sensitivity strip. 

In the case of the dipole magnetic field, the differences be- 
tween the magnetic and nonmagnetic rb(rs)-relations become 
appreciable only in a very strong magnetic field, B > 10^* G, 
and if the internal temperature is relatively low, Tb < 3 x 10^ K. 
In this case, Ts is noticeably increased by the magnetic field. 
The increase of a mean effective temperature is however quite 
moderate in a comparison to the effect of such a strong field on 
the local temperature value: a significant increase of at the 
magnetic pole and a sharp decrease at the equato r (see, e.g. 



I l200lt IPotekhinetaLl l2003V For lower 
magnetic field and higher temperature, the difference in the 
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Fig. 2. The dependence of the mean effective temperature on 
B for the dipole (dot-dash lines) and small-scale (solid lines) 
fields. The curves are shown for the accreted and iron envelopes 
and the internal temperature T\, = 10^ K and 10** K (marked 
near the curves). 

magnetic geometry does not yield a considerable departure in 
the rb(rs)-relation. 

The difference in a magnetic configuration can manifest it- 
self at B > lO"' G. The small-scale field, as a rule, results 
in a lower at a given B compared to the dipole field. The 
difference in for the small-scale and dipole fields depends 
generally on Th, B, and the mass of the accreted light-element 
material (see Fig.|2j, but never exceeds 20%. However, this dif- 
ference in Ts can result in a factor of ~ 2 in the luminosity. In 
all considered models, the uncertainty in a chemical compo- 
sition of the envelope (light versus heavy elements) causes a 
much larger difi'erence in than the uncertainty in the mag- 
netic field geometry (dipole versus small-scale). 

To estimate the effect of a small-scale magnetic field on the 
rb(rs)-relation, one can use the fit ting formulae derived for the 
dipole field (see the Appendix of IPotekhin et al.lf2003l) where 
cos^fl and sin^ should be replaced by 1/3 and 2/3, respec- 
tively. Such fit reproduces the present numerical results with 
an accuracy < 10%. 

4. Discussion 

We have considered the relation between the surface and inter- 
nal temperature in neutron stars in the case where small scale 
magnetic fields near the surface are stronger than a large scale 
(e.g., dipole) field. Calculations show that the difference in the 
Th(T s)-Telation between the stars with small-scale and dipole 
fields is not very significant, although, generally, it should be 
taken into account if highly accurate thermal luminosity calcu- 
lations are required. This can be important for the interpretation 
of future high-precision measurements of the neutron star ther- 
mal radiation. 



For both the small-scale and dipole fields, the departure 
from the thermal structure of a nonmagnetic neutron star is rel- 
atively small, if the typical field strength B < 10'"^ G. This is 
caused by the fact that the increase of thermal insulation near 
the region where the field lines are tangential is well compen- 
sated by the decrease of the insulation in the region where the 
field is normal to the surface. As seen from Fig.|2] for the dipole 
field this compensation is generally more efficient than for the 
small-scale field. 

For a stronger field, an increase of the thermal conductivity 
along the field lines due to the magnetic quantization effects 
turns out to be so strong that it cannot be fully compensated by 
a decrease of the conductivity in the perpendicular direction. 

Our results indicate that the geometry and strength of the 
magnetic field are likely unimportant for the average ther- 
mal structure of neutron stars, if the field strength is moder- 
ate (B < 10'"^). Therefore, the magnetic field of a moderate 
strength can affect the neutron star cooling via Joule heating 
rather than via the rb(rs)-relation. Joule heating is of particu- 
lar importance at the late evolutionary stage and can maintain 
a relatively high surface temperature > 10^ K during a very 
long time comparable to the decay time of the magnetic field 
(Miralles etal., 1998). 
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